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Abstract 

Recently it was established that a certain integrable lonq-ranqe spin chain describes the 
dilatation operator of A/" = 4 gange theory in the su(2) sector to at least three-loop order, 
while exhibiting BMN scaling to all orders in pertnrbation theory. Here we identify this 
spin chain as an approximation to an integrable short-ranged model of strongly correlated 
electrons: The Hnbbard model. 



1 Introduction 


Recently it was discovered that the planar one-loop dilatation operator of snpersymmetric 
A/” = 4 gange theory is completely integrable pj |2j. This means that its spectrum may 
be exactly determined in the form of a set of non-linear Bethe equations. Evidence was 
found that this integrability is preserved beyond the one-loop approximation, and it was 
conjectured that the dilatation operator might be integrable to all orders in perturbation 
theory [3]. Given the usually benign, analytic nature of planar perturbation theory, one 
may then even hope for the theory’s complete large N integrability at all values of the 
Yang-Mills coupling constant. 

Deriving the dilatation operator from the held theory, and subsequently demonstrat¬ 
ing its integrability, is not easy. The three-loop planar dilatation operator in the max¬ 
imally compact su(2|3) sector was found by Beisert, up to two unknown constants, by 
algebraic means in j^. These constants could later be unequivocally hxed from the 
results of a solid held theory calculation of Eden, Jarczak and Sokatchev [H]. This ba¬ 
sically completely determines the planar dilatation operator in this large sector up to 
three loops. Its restriction to su(2) agrees with the original conjecture of j2]. Three-loop 
integrability in su(2) was then demonstrated in [H] by embedding the dilatation opera¬ 
tor into an integrable long-range spin chain due to Inozemtsev, and a three-loop Bethe 
ansatz was derived. 

The Inozemtsev spin chain exhibited a four-loop breakdown of BMN scaling jHj. This 
scaling behavior seemed, and still seems, to be a desirable, albeit unproven, property of 
perturbative gauge theory. Mainly for that reason an alternative long-range spin chain, 
dihering from the Inozemtsev model at and beyond four loops, was conjectured to exist 
in [Zj. Its construction principles were an extension of the ones already laid out in [Hj: 
(1) Structural consistency with general features of Yang-Mills perturbation theory, (2) 
perturbative integrability and (3) qualitative BMN scaling. The model’s Hamiltonian is 
only known up to hve loops, and increases exponentially in complexity with the loop or¬ 
der. In striking contrast, a very compact Bethe ansatz may be conjectured for the model 
and shown to diagonalize the Hamiltonian to the known, hfth, order. The conjecture 
reads 
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where the rapidities Uk = u{pk) are related to the momenta pk through the expression 


u{Pk) = ^ cot Y + 8^2 gin2 ^ ^ ^2) 

and the energy should be given by 



This Bethe ansatz should yield the anomalous dimensions A of su(2) operators of the 
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form 


+ ... , where A{g) = L + E{g) with / = ^ = A. (4) 

The dots indicate all possible orderings of the partons Z and X inside the trace. This 
mixing problem is diagonalized by the spin chain Hamiltonian, where we interpret Z 
as an up-spin | and X as a down spin j. L is the length of the spin chain, and M 
the number of magnons J,. These are the elementary excitations on the ferromagnetic 
vacuum state | || ... |t) which should be identified with the gauge theory’s BPS state 
Tr Z^. To leading one-loop order the spin chain Hamiltonian coincides with the famous 
isotropic nearest-neighbor Heisenberg XXX spin chain P, and the corresponding Bethe 
ansatz is obtained by taking the g ^ 0 limit of o,©,©. See also P for a detailed 
explanation of the long-range spin chain approach to gauge theory. 

The higher-loop Bethe ansatz ©,©,© has many intriguing properties [7j, and it is 
suspicious that it should not have already appeared before in condensed matter theory. 
It is equally curious that the Hamiltonian should be so complicated, see [Zj, to the point 
that it is unknown how to write it down in closed form. Finally, and most importantly, 
the Bethe ansatz is expected to break down at wrapping order, i.e. it is not believed to 
yield the correct anomalous dimensions A starting from 0{g‘^^). This suggests that the 
asymptotic Bethe ansatz 0,0,® is actually not fully self-consistent at hnite L and 
9 ^ 0 - 

Nearly all work in solid state theory on the Heisenberg magnet has focused on the 
antiferromagnetic vacuum and its “physical” elementary excitations, the spinons. The 
only notable exceptions seem to be two articles of Sutherland and of Dhar and Shastry 
[Tn] , where it was noticed that the dynamics of magnons in the ferromagnetic vacuum is 
far from trivial. This was later independently rediscovered and extended in the A/” = 4 
context in mnn. In gauge theory the BPS vacuum is very natural, but it should be 
stressed that all states are important. In particular, it is interesting to ask what is the 
state of highest possible anomalous dimension. This is precisely the antiferromagnetic 
vacuum state, where M = L/2 and E{g) in © should be maximized. Contrary to the 
BPS state | || ... ||) this state is highly nontrivial, as the Neel state | ti-Ti- • • • Ti-) is 
not an eigenstate of the Heisenberg Hamiltonian. This problem was solved for = 0 in 
the thermodynamic limit L —00 in 1938 by Hulthen ^21 using Bethe’s ansatz. 

Like the BPS state, the antiferromagnetic vacuum state is of very high symmetry. 
It should therefore also be of great interest in gauge theory. Let us then use the BBS 
Bethe ansatz 0,0,© and compute the higher-loop corrections to Hulthen’s solution. 
As the computation is done in the thermodynamic limit the BBS equations are perfectly 
reliable. The one-loop solution may be found in many textbooks. It is particularly 
well described in the lectures jHj. Adapting it to the deformed BBS case is completely 
straightforward. We will therefore mostly skip the derivation, referring to for details, 
and immediately state the result for the energy of the antiferromagnetic vacuum: 

E{0)=L f du ^ E{g) = L f du p{u) (, (5) 

i_oo U^ + J J-oc \x+{u) X {u)J 
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where the auxiliary spectral parameter a; [Zj is given by 


x{u) 


u 
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with x^{u) = x{u ± -). 


( 6 ) 


Here p{u) is the thermodynamic density of (magnon) excitations. It is found from solving 
the Bethe equations, which turn at L —> cx) into a single non-singular integral equation 
for p{u)\ 


dp{u) 

du 


2n p{u) -|- 2 
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where the derivative of the momentum density is, with u± = u ± |, given by 


dp{u) 

du 
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We notice that the r.h.s. of Hulthen’s equation m does not depend explicitly on the 
coupling constant g (since the S-matrix of the BDS Bethe equation, i.e. the r.h.s. of 
m does not look different, in the n-variables, from the one of the Heisenberg model). 
Furthermore, the kernel of the integral equation is of difference form and the integration 
range is inhnite. The equation may therefore immediately solved for p{u), for all g, by 
Fourier transform: 


p{u) 


9=0 


2 cosh Tiu 


, , dt cos (tu) Jo(\/2gt) 

PW = TtZ - 


27r 


cosh (I) 

Plugging this result into the energy expression © one hnds 

4 p dt Jo{V2gt) Ji{V2gt) 


(9) 


E{0) = L2 log2 


E{g) = L 
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where Jo{t), Ji{t) are standard Bessel functions. 

Now, it so turns out that the expressions for p{u) in 0 and E{g) in are very 
famous results in the history of condensed matter theory. The latter is, up to an overall 
minus sign, identical to the ground state energy of the one-dimensional Hubbard model 
at half hlling. It was shown to be integrable and solved by Bethe Ansatz in 1968 by Lieb 
and Wu ini. Since then a very large literature on the subject has developed. For some 
good re- and overviews, see ng. na . The Hubbard model is not quite a spin chain, 
but rather a model of Nq itinerant electrons on a lattice of length L. The electrons are 
spin-| particles. Due to Pauli’s principle the possible states at a given lattice site are 
thus four-fold: (1) no electron, (2) one spin-up electron |, (3) one spin-down electron |, 
(4) two electrons of opposite spin Hubbard’s Hamiltonian reads, in one dimension 
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The operators c\^ and are canonical Fermi operators satisfying the anticommutation 
relations 


{4,cp}=0, (12) 

^ij ^UT ■ 

We see that the Hamiltonian consists of two terms, a kinetic nearest-neighbor hopping 
term with strength t, and an ultralocal interaction potential with coupling constant U. 
Depending on the sign of U, it leads to on-site attraction or repulsion if two electrons 
occupy the same site. 

Comparing the BDS result m with the result of Lieb and Wu for the ground state 
energy of the half-hlled band, where the number of electrons equals the number of lattice 
sites, i.e. Nq = L, we see that the two energies coincide exactly under the identihcation 


{Ci,0-5 

c] T_} = 


1 



(13) 


This leads us to the conjecture that the BDS long-range spin chain, where, by construc¬ 
tion, g is assumed to be small, is nothing but the strong coupling limit of the Hubbard 
model under the identihcation (HI. In the following we will show that this is indeed the 
case, even away from the antiferromagnetic ground state. In fact, we shall demonstrate 
that it is exactly true at finite L up to 0{g‘^^) where the BDS long-range chain looses its 
meaning. This will, however, require the resolution of certain subtleties concerning the 
boundary conditions of the Hamiltonian (ED. As it stands, it will only properly diago¬ 
nalize the BDS chain if the length L is odd. It the length is even, we have to subject the 
fermions to an Aharonov-Bohm type magnetic hux f. The Hamiltonian in the presence 
of this flux remains integrable and reads 


H = 


'^g i=l a=T,i 
where the twist is given by^ 


Lj 1 ^ 

i=l a=T,i ^ i=l 


J 




0(T f 1 cr T) i ) 

0 = 0 for L = odd and 


0 = — for L = even. 
^ 2L 


( 14 ) 


( 16 ) 


An alternative way to introduce the Aharonov-Bohm flux is to perform a suitable gauge 
transformation and to thereby concentrate the magnetic potential on a single link, say 
the one connecting the L’th and the hrst site. It is then clear that considering a non-zero 
flux amounts to considering twisted boundary conditions for the fermions. 

The vacuum of the Hamiltonian m is the empty lattice of length L. Here the 
elementary excitations are up (t) and down (j) spins. Two electrons per site (|) are 

^For odd L the twist could alternatively be chosen as any integer multiple of while for even L 
any odd-integer multiple of ^ is possible. A compact notation which does not distinguish the cases L 
odd or even is 0 = ■ 
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considered a bound state of elementary excitations. These constituents of the bound 
states are repulsive (as > 0). For our purposes it is perhaps more natural to consider 
the BPS vacuum: 


\Z^) = I II .. 41) = |0) (16) 

We may then perform a particle-hole transformation on the up-spin electrons. 
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(17) 

(18) 


Now single up-spins (|) are considered to be empty sites, while the elementary excitations 
are holes (o) and two electrons states (4- In the condensed matter literature, such a 
transformation is often called a Shiba transformation and it is known to reverse the sign 
of the interaction. The standard Shiba transformation contains an alternating sign in the 
dehnition of the new creation/annihilation operators, designed to recover the hopping 
term, at least in the periodic case. The price to pay is that for odd lengths the sign of the 
hole hopping term will change on the link connecting the last (L’th) and the hrst site. 
In other words, the particle/hole transformation introduces an extra flux of vrL seen by 
holes. Since we prefer to distribute this twist uniformly along the chain, we remove the 
signs in the dehnition of the hole operators^ and put 

Ci,o = , (19) 

Ci,i = Ci^i , d,i = • (20) 

Under the particle/hole transformation, the charge changes sign and the corresponding 
hopping terms get complex conjugated. An extra minus sign comes from the reordering 
of the hole operators. Therefore we may write the Hamiltonian in its dual form 


H = 


L 
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'' ^ 1=1 cr=o,l 






Ed 


4 c 




, 1^81 


( 21 ) 


2=1 


where while 0o = tt — (j)^. Comparing the two expressions and m we 

conclude that under the duality transformation, the Hamiltonian (HI transforms as 

^^ ( 22 ) 

As predicted, the sign of the interaction changes upon dualization. The effect is that 
holes o and states with two electrons per site | attract each other and form bound states 
I, the magnons. 

^This amounts to a gauge transformation. 
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2 Effective Three-Loop Spin Hamiltonian 

In this chapter we will explicitly demonstrate that the Hubbard Hamiltonian m gen¬ 
erates at small g the three-loop dilatation operator of A/” = 4 gauge theory in the su(2) 
sector jS]. The BDS long-range spin Hamiltonian j2j is thus seen to emerge as an effec¬ 
tive Hamiltonian from the underlying short-range system. Note that the small g limit, 
relevant to perturbative gauge theory, corresponds, via m, to the strong coupling limit 
U ^ oo oi the Hubbard model in condensed matter parlance. 

Our claim may be verihed immediately to two-loop order, using well-known results 
in the literature. Klein and Seitz HE] proposed the strong-coupling expansion of the 
half-hlled Hubbard model to 0{g'^). The two-loop result 0{g'^) was later conhrmed by 
Takahashi Hg. In fact, eq. (2.15) of his paper^ precisely agrees with the two-loop piece of 
the BDS Hamiltonian (and therefore with two-loop gauge theory [3]) under the parameter 
identihcation OBJ. Eq. (2.15) of pH] also contains certain four-spin terms which only 
couple, since our system is one-dimensional, to a length L = 4 ring. These are a first 
manifestation of certain unwanted terms which we need to eliminate by appropriate 
boundary conditions and twisting, see (Cl,®, to be discussed in more detail below. 

When one now turns to the three-loop 0{g^) result of Klein and Seitz as ob¬ 
tained in 1972, one unfortunately finds that their effective Hamiltonian disagrees with 
the BDS Hamiltonian at this order. We have been unable to hnd a later paper in the 
vast condensed matter literature on the subject which conhrms or corrects their 33 year 
old calculation. We have therefore decided to check their computation in detail. And 
indeed we found a mistake, see below. Correcting it, we reproduce the planar three-loop 
su(2) dilatation operator HUE], see (j2H|) , (d below. 

For the remainder of this chapter it is convenient to use dH and rewrite it (see 
Appendix IA.3I for a discussion about the relevance of the twist factors in computing the 
effective Hamiltonian) in the form: 

L L 

T^Hubbard ^ ^ T Xji) -\- tU ^ ^ ■ (^2) 

i<j i=l 

where Xij = ^ • 

2.1 Generalities 

The Hamiltonian (d consists of two parts: A hopping term involving the coefficients 
tij, and the atomic part. The latter is diagonalized by eigenstates describing localized 
electrons at sites Xj. The ground-state subspace £q of the atomic part is spanned by 
4 ti 4 t 2 • • • TL-i |0)- Here we are interested in the limit of large H, with t staying 
relatively small. The atomic part tends to localize the electrons, while some hopping 
may still occur. At low temperatures this corresponds to small fluctuations around To 
states, since each hopping of the electron from one site to another is suppressed by a 

^Incidentally, this is the famous paper where the next-nearest neighbor correlation function of the 
Heisenberg antiferromagnet was first obtained. We took this as a hint that the half-filled Hubbard model 
“knows” something about long-range deformations of the Heisenberg model. 
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factor of order of 1/U. One can now pose the question whether it is possible, for large 
U and low temperature, to hnd an effective operator h acting in Sq whose eigenvalues 

h|0) = E|0) (24) 

are the same as for the one of the Hamiltonian (ED: 

H\^l^) = E\^P). (25) 

The answer is to the positive and has a long history [20] . A formal and rigorous treatment 
of this subject is presented in appendix lAl 

It is however instructive to discuss (I24|l in a more heuristic way. It is obvious that 
the effective Hamiltonian h must properly include the hopping effects. On the other 
hand it acts only in a subspace of the full state space, where conhgurations with double 
occupancies are projected out. This means that (1^ should describe processes with 
virtual intermediate states, corresponding to electrons hopping from site i to site j and 
subsequently hopping back. Since every nearest-neighbor hopping is suppressed by 1/f/ 
it is clear that the 1/U expansion of h will result in increasingly long-range interactions. 
What kind of terms may appear in hi A hrst guess leads to products of hopping operators 
Xij with the condition that they will not move states out of the space £q. Since X^j 
annihilates an electron at site j and creates a new one at i, we see that only such products 
of Xij operators are allowed which result in the same number of creation and annihilation 
operators at a given lattice site. Since each product of creation and annihilation operators 
may be represented in terms of 5u(2) spin operators, we conclude that the effective 
Hamiltonian (1^ must be of spin-chain form! 

2.2 Three-Loop Result 

We have used perturbation theory for degenerate systems (see appendix na where also 
some details of the computation scheme are explained) to derive the effective Hamiltonian 
to three loops. The result up to sixth order (i.e. three loops) for the formal perturbation 
theory expansion is found in in appendix IA.2I It may be shown to be completely 
equivalent to the expansion obtained by Klein and Seitz in El- 

The formal expansion is then converted into a diagrammatic expansion, see again 
appendix IA.21 We agree with Klein and Seitz with all perturbation theory diagrams 
up to sixth order as presented in their paper, except that we hnd that they missed a 
few diagrams (of type ’a’ as in Fig 5. of their paper). These are the following diagrams 
(summation over i is understood): 




+ 


i i+1 i+2 i+3 
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where ... means arrow-reversed diagrams. 

We have conhrmed all diagram evaluations performed in except for the contri¬ 
bution of the diagrams of type / in equation (C3) of the mentioned paper, where there 
is an overall factor of 16 missing. We believe this to be a typographical error. There is 
however also an additional contribution from the mentioned four diagrams which were 
not included in their computations. Explicit calculation shows, that the missing terms 
yield 

6 

Ut{lQAi - 4 A 2 + 2 B 3 - 2 Bi - 2 B 2 ), (26) 

where 

~ 'Pi,i+s) ! -^1 = ~ 'Pi,i+l'Pi+2,i+3) , 

-®2 = “ 'Pi,i+2'Pi+l,i+3) , B3 = ~ 'Pi^i+3'Pi+l^i+2) , 



and P is a spin permutation operator. Correcting the result of Klein and Seitz we hnd 


h 



tUAi + 
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Upon putting U 
form 


with 


VI + — _K 

9 ’ ^ “ V 2 g 
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tU{B 2 - B 3 ) 


(27) 

and after some simple algebra one rewrites (EID in the 


L 

h = ^^(^2 + + ...), 

i=l 


(28) 


^2 — 2 ~ 'TCi+l) 5 

^4 = —(1 — O’ifi+l) + —(1 — Cri(Tj+2) , 

15 3 1 

he = —(1 — —-(1 — + ^(1 ~ di <? i - i - 3 ) 

~ ~ ^i+l^i+ 2 ) 

o 

+ -(1 ~ <Taj+2)(l — Uj+iUj+s). (29) 

O 


This is indeed the correct planar three-loop dilatation operator in the su(2) sector of 
M = 4 gauge theory (3]. It is fascinating to see its emergence from an important and 
well-studied integrable model of condensed matter theory. 






3 Lieb-Wu Equations 

The Hamiltonian (ED was shown to be integrable and diagonalized by coordinate Bethe 
ansatz in • For a pedagogical treatment see . This required finding the dispersion 
relation of the elementary excitations | and J, and working out their two-body S-matrix. 
It is indeed a matrix since there are two types of excitations, hence their ordering matters. 
The scattering of two up- or two down-spins is absent, as identical fermions behave like 
free particles. The scattering of different types of fermions is non-trivial due to their 
on-site interaction. After working out the S-matrix one needs to diagonalize the multi¬ 
particle system by a nested Bethe ansatz. The result of this procedure, generalized to 
the case with magnetic flux, yields the Lieb-Wu equations: 


^^qr^L _ TT Uj - V2gsm{qn + (l))-i/2 _ 

I I ^ ^ ^ 

j=i “i - V2gsm{qn + 0) + i/2 

A Uk - V2gsm{qn -f- 0) + i/2 _ A Uk -Uj+i 
Uk - V2g sm{qn -f 0) - i/2 Uk-Uj-i' 


(30) 

(31) 


where the twist is given"^ in ED and the energy is 


E 


A 
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L 

cos(g„ -F 0) . 

n=l 


(32) 


Here we have already specialized to the half-filled case with Nq = L fermions and M < 
Ll2 down-spin fermions (there are thus L — M up-spin fermions in the system). 

This form of the Hubbard model’s Bethe equations if very convenient for demonstrat¬ 
ing rather quickly that the 0 limit yields the spectrum of the Heisenberg magnet. 
In fact, the Lieb-Wu equations decouple at leading order and become 




ivi . , 

n uj - V2 

Uj +i/2 


(^1 -h 0(5')) , 


/ Uk -F i/2 \^ _ A Uk -Uj + i 
\uk-i/2) ^}^Uk-Uj-i' 


n = 1,..., L 
k = 1,...,M 


(33) 

(34) 


Eqs. dni are already identical to the ones of the Heisenberg magnet (see e.g. PiEl)- 
The r.h.s. of is, to leading order 0(5°), the eigenvalue of the shift operator of the 
chain (again, [I],[I1]). In gauge theory we project onto cyclic states, so we may take the 
eigenvalue to be one, and solve immediately for the L momenta qn to leading order: 

- r 2TI 

^ ^ _ 1 ) + n = l,...,L. (35) 

'*^The Lieb-Wu equations for arbitrary twist are given in appendix |0 
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But now we have to find the energy. Plugging the result into the expression (|S2D 
conveniently eliminates the spurious 0{l/g) term in the energy. We therefore need to 
find the 0{g) corrections to the momenta in ()d5|l from ()d()|l . Luckily, this is a linear 
problem; solving it one computes the 0{g^) term of as 

M 

^ = E lUTT + 0(9) ■ (36) 

k=l “fc ^ 4 

which is the correct expression for the energy of the Heisenberg magnet. 

The starting point for the small g expansion of the Lieb-Wu equations are therefore 
Bethe’s original equations (ini,® in conjunction with the free particle momentum 
condition (jSSD- It is interesting that all non-linearities are residing in the one-loop Bethe 
equations (j^H). Once these are solved for a given state, the perturbative expansion is 
obtained from a linear, recursive procedure. It allows for efficient and fast numerical 
computation of the loop corrections to any state once the one-loop solution is known. 
A simple tool for doing this with e.g. Mathematica may be found in Appendix El along 
with a similar tool for the perturbative evaluation of the BDS equations. 

We have applied this perturbative procedure to alB (cyclic) states of the BDS chain 
as recorded, up to five loops, in Table 1, p.30 of [7j. The (twisted) Lieb-Wu equations 
dsni),®,® perfectly reproduce the energies of this table. 

We found that that our version of the Hubbard model precisely agrees in all investi¬ 
gated cases with the results of the BDS ansatz up to and including the (L — l)-th loop 
order. On the other hand, invariably, at and beyond L’th order of perturbation theory 
(corresponding to the 0{g‘^^~^) terms in the energy E) the predictions of the two ansatze 
differ. See chapter E for some concrete examples. 

It is also interesting to record the effects of the twists on the perturbative spectrum. 
A first guess might be that they should only influence the spectrum at and beyond 
wrapping order 0{g‘^^~‘^), when the order of the effective interactions reaches the size of 
the ring, and the system should become sensitive to the boundary conditions. In actual 
fact, however, one finds that the twists generically influence the spectrum starting at 
already 0{g^~‘^). This is the phenomenon of demi-wrappings. The Hubbard model 
at small g behaves effectively as a long-range spin chain due to the virtual “off-shell” 
decomposition of the magnon bound states | (which are sites occupied by a down-spin 
but no up-spin) into holes o (empty sites) and double-occupied sites |. The power of 
the coupling constant g counts the number of steps a hole o or |-particle is exercising 
during its virtual excursion, see also the discussion in section 12.11 We now observe 
that starting from at 0{g^~^) the excitations o,^ can (virtually) travel around the ring, 
and the amplitudes start to depend on the boundary conditions! A similar distinction 
between wrappings and demi-wrappings was qualitatively discussed in a recent paper on 
this subject |22]- Our procedure of twisting eliminates the demi-wrappings. Interestingly, 
this seems to leave no further freedom at and beyond wrapping order, at least in the 
context of our current construction. 

The Lieb-Wu equations in the form (isni),®,® are very useful for the analysis of 
chains of small length. They are far less convenient in or near thermodynamic situations, 

®The only exception are certain singular three-magnon states which require a special treatment. 
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i.e. when L —>■ cx). The reason is the large number of momenta one has to deal with. 
In (|?H) we have written a dual form of the Hamiltonian (fT^ . Accordingly, we may write 
down the corresponding set of dual Lieb-Wu equations: 



hi “li “ '/2gsin{q„ -<!>) + i/2 ’ 



(37) 



M/j Uj I 

j=l J 


k = l,...,M (38) 


where the energy is now given by 



(39) 


Again, we have specialized to the case of half-hlling. A particular feature of the dual 
Hamiltonian m is that the twist is different for the two components. We are therefore 
led to use the Lieb-Wu equations for generic twist which are written down in Appendix 
O This explains the minus sign in the right hand of = —1. 

Note that cj) —^ —cj) is a symmetry of the equations (but not of the solutions), as we may 
change u —> —u and q —>■ —q. Note also that therefore the set of L -|- 2M momenta 
idn, —qn) corresponds to the L + 2M solutions of the first Lieb-Wu equation (jHIIl)- 

4 Magnons from Fermions 

In chapter El we proved, to three-loop order, that the Hamiltonian of the BDS long- 
range spin chain emerges at weak coupling g from the twisted Hubbard Hamiltonian as 
an effective theory. Pushing this proof to higher orders would be possible but rather 
tedious. Note, however, that the BDS Hamiltonian is, at any rate, only known to hve- 
loop order [Zj. What we are really interested in is whether the Bethe ansatz (HD,©,®, 
which was conjectured in [2] , may be derived from the Bethe equations of the Hubbard 
model, i.e. from the Lieb-Wu equations of the previous chapter. We will now show that 
this is indeed the case. The derivation will first focus on a single magnon fsection l4.1|l . 
where it will be shown that the magnons | of the long-range spin chain are bound states 
of holes o and double-occupations as is already suggested by the perturbative picture 
of chapter El It will culminate in EM where we demonstrate that the bound states alias 
magnons indeed scatter according to the r.h.s. of (P). An alternative proof may be found 
in appendix lEl 

Unlike the BDS long-range spin chain, the twisted Hubbard model is well-defined 
away from weak coupling, and actually for arbitrary values of g. An important question 
is whether the twisted Hubbard model allows to explain the vexing discrepancies between 
gauge and string theory 123 E]. Unfortunately this does not seem to be the case. We 
have carefully studied the spectrum of two magnons in section 14.21 and find that their 
is no order of limits problem as the coupling g and the length L tend to infinity while 
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gjL stays finite. The scattering phase shift indeed always equals the one predicted by 
the BDS chain. 

The Hubbard model contains also many states which are separated at weak coupling 
by a large negative energy gap 0{—l/g‘^) from the magnons. This may be seen from the 
expression (jsni). For solutions with real momenta Qn the cosine is bounded in magnitude 
by one, and the constant part —M/g"^ cannot be compensated. These states are composed 
of, or contain, holes o and double-occupations \ which are unconhned, i.e. which do not 
form bound states. Their meaning will need to be understood if it turns out that the 
A/" = 4 gauge theory’s dilatation operator can indeed be described by a Hubbard model 
beyond the perturbative three-loop approximation. In fact, it is clear from the expression 
for the anomalous dimension A = L + g'^ E{g) in (Q that each unconfined pair (o, |) 
shifts the classical dimension and thus the length down by one: L ^ L — 1. Is this a 
hrst hint that the perturbative su(2) sector of A/” = 4 gauge theory does not stay closed 
at strong coupling, as was argued in [21]? 

4.1 One-Magnon Problem 

Let us then begin by studying the case of M = 1 down spin and L — 1 up spins, see 
[IB],|22j. Clearly it is easiest to use the dual form of the Lieb-Wu equations (HII),®,®. 
In the weakly coupled spin chain we have only L states, while in the Hubbard model 
we have states. This is because one down spin J, is composed of one hole o and one 
double-occupation |. If we project to cyclic states, as in gauge theory, only one of the 
L states survives, namely the zero-energy BPS state. However, in order to derive the 
magnon dispersion law, we will not employ the projection for the moment. This way the 
magnon can carry non-zero momentum and energy. In the Hubbard model the magnon 
should be a o— | bound state, and we therefore make the ansatz (with (3 > 0 and g > 0): 

gi-0=| + g-l-i^, q^-^=^ + q-il3. (40) 

Here qi and q 2 are the quasimomenta of the o and the | particles. They are complex, 
where the imaginary part (3 describes the binding. Adding the real parts gives the 
momentum 2q of the magnon. The dual Lieb-Wu equations for one magnon, where we 
only have a single rapidity u, read 

gigiL ^ u-^/2gsin{qi -^)-i/2 ^ u - V 2 gsm{q 2 - (j)) -f/2 

u — \/2gsm{qi — 0 ) -|- z/2 ’ u — v^gsin(g 2 — 0) + ' 

u — ^/2g sin(gi — 0) -|- i/2 u — \/ 2 gsm.{q 2 — (f)) -\- i/2 ^ 

u — \/2gsm.{qi — </>)— i/2 u — V2g sin(g 2 — (j)) — i/2 

By multiplying, respectively, the left and right sides of the two equations in dn} 
using (HOD,® we derive 

^ q=ln (u = 0,1,...,L-1). 

Jj 

This is just the statement that the magnon is free (there is nothing to scatter from) and 
its momentum p := 2q is quantized on the ring of length L. Furthermore we can rewrite 


( 41 ) 

( 42 ) 
and 

( 43 ) 
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(PT|l as 


V2g sin (gi ,2 - 0) - m = - cot 


qi,2L 

2 


(44) 


Decomposing into real and imaginary parts we find, using the twist Cl, 


sinh (/?) 


1 1 
2V2g sin (g) 


tanh (/? L) , 


(45) 


and® 

u = V2g CO. (g) co.h [0) + ■ (46) 

By analyzing iSl) we may now discuss the existence of bound states. We see that for 
large L, where tanh {f3L) —1, we have, for given mode number n, exactly one^ solution 
with /3 > 0 for all values of g > 0. We also see that there is only one way to take the 
thermodynamic limit, independent of g: 


sinh (/3) ~ 


1 1 
2 \/ 27 t ng/L ' 


(47) 


But this means that there is also only one way to take the BMN scaling limit, where 
g,L —>■ oo with g/L kept finite. 

Let us then work out the energy of the magnon with momentum p = 2q. The 
exponential terms tanh(/3L) ~ 1—may clearly be neglected at large L for arbitrary 
values of g, and we immediately find the dispersion law 


E = -4 + ^^ sin cosh(/3) = + %g‘^ , (48) 

r 9 ^2/ g^y 2 

which is exactly the BDS result ©! Likewise, again dropping the exponential terms 
from the rapidity relation (gni), we find the BDS result 0 for the dependence of the 
rapidity u{p) on the momentum Pk = P = 2g. 

Note that our derivation only assumed the thermodynamic limit; it did not assume 
weak coupling. If the coupling g is weak we may in addition deduce from (gSl) that the 
binding amplitude (3 diverges logarithmically as ~ — logg. We may then deduce that 
the exponential terms we dropped are 


e-2/3L 




(49) 


and therefore should be interpreted as 0{g‘^^) wrapping corrections. 

We just showed that L of the states of the Hubbard model’s M = 1 states can be 
interpreted as magnons. The remaining L{L — 1) states should correspond to solutions 

®The sign of the second term in iHHIl may be changed by choosing a different gauge for the twist. 
This type of gauge dependence should not appear in physical observables such as the energy. 

^Actually, if g becomes of the order of L such that gjL is larger than a certain threshold value, the 
bound state is lost. An additional real solution, c.f. appendix IHI will appear. 
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where the momenta gi, q 2 are real, i.e. these are not bound states. Among these, L — 1 
states are cyclic. The unbound states are found as follows. We make the ansatz 

71 71 

gi-0 = - + g + &, q 2 -(j) = - + q-b, (50) 

which is completely general except for the assumption that q and b are real. The twisted 
dual Lieb-Wu equations (gH),® still apply, and, using the same multiplication trick as 
before we hnd again (HHl) . This is merely the statement that that the total momentum 
is quantized on the ring of length L. The Lieb-Wu equations m now read 


V29 cos ((/ ± 6) - ti = 1 cot !((/ ± 6 + ^ + ^) A 


(51) 


Let us hrst consider the unconhned cyclic states, i.e. the case of mode number n = 0 
in (gni), hence q = 0. We can then immediately hnd the energy of such states to-be from 

da 

B = -4. (52) 

which is seen to not depend on b. But can we really hnd values for b which satisfy the 
Lieb-Wu equations (ED? How many solutions of this type do we have? The answer is 
easily found from subtracting either side of the two equations in (ED- This yields the 
consistency condition 

lcot((|+.#. + (,)^) = loot ((! + .#.-(,) I) . (53) 


Now it is very easy to show that there are precisely L — 1 solutions of this equation: 

TT 

b = —m with m = l,...,L —1. (54) 

Jj 

Therefore, the M = 1 cyclic unconhned {L — 1-fold degenerate) states resemble zero- 
energy “BPS states” with exact scaling dimension A = L — 1. However, see appendix 

o 

Finally, let us study the number of states and the dispersion law of the unconhned 
states carrying non-zero total momentum p = 2q = ^n, cf. (ED- ^^at by 

eliminating u from ED that 

sin (6) = —^-- tan(6L) , (55) 

2^/2g sm (g) ' ' 

which turns out to just be the analytic continuation of dUD. It is not hard to prove that 
there are indeed generically L — 1 solutions for each value of the L — 1 non-zero values of 
g. This yields (L — 1)^ states. Therefore, adding these to the L — 1 cyclic real solutions, 
and the L bound states, we have accounted for all of the states of the M = 1 problem. 
In appendix Owe investigate the energy of the real solutions in the large g and large L 
limit. In the limit L —>■ 00 , the solutions of (ED become dense on the interval (0, 27r), so 
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for any value of the magnon momentum p = 2q we have a continuum of states, whose 
energies vary continuously. It is not clear how one would interpret these states in the 
context of the gauge theory, or, more generally, the AdS/CFT correspondence. It is 
possible that we need a model encompassing all the sectors of the gauge theory to be 
able to draw some conclusion about the large g limit. 

Let us now turn to the mutual scattering of our magnons; hrst for two, and then for 
arbitrarily many. We shall hnd that the scattering is, up to exponential terms, indeed 
given by the r.h.s. of (HD- 


4.2 Two-Magnon Problem 

The result of the previous section does not bode well for the hope expressed in the last 
chapter of [7j that wrapping might explain the discrepancies between gauge and string 
theory. This would require an order of limits problem as one takes the coupling g and 
the length L large. It is certainly not seen on the level of bound state formation, recall 
(USD- However, one might still hope that the magnons constructed in the last section 
might somehow scatter in distinct ways at weak and strong coupling. By considering the 
M = 2 two-magnon problem we will now show that, unfortunately, this is not the case. 
It therefore seems that the AFS string Bethe ansatz [21] cannot be obtained from the 
twisted Hubbard model, at least not in the current version. 

We proceed much as before, making the appropriate ansatz for two holes and two 
double-occupancies (with /3 > 0 and g > 0) bound into two magnons with momenta 
p = 2q and —p = —2q: 

71 71 

gi-(/> = - + q + ip, q3-(j)=^ + q-ip, (56) 

g2-0 = + q4-(j) = -'^-q-i(^- 

We derive (for simplicity assume L = 1 mod 4, which allows to assume that the two 
rapidities obey ui = —U 2 ) from the dual Lieb-Wu equations 


sinh (/3) 


cosh (/3) cot (g) sinh {(3 L) 

2^/2g cos (g) cosh (/3) (cosh (/3 L) -f sin (gT)) — cos (gT) 


(57) 


While looking superhcially different, this agrees precisely, up to exponential terms, in 
the L —> cx) limit with (jUj). It is crucial to note that, as in the previous M = 1 case, 
there is only one way to take the BMN scaling limit, which yields again dlZD. Likewise, 
we hnd, up to exponential corrections, for the rapidity 

= ^ + g‘^ + g‘^ cos (2g) cosh [2(3) = ~^ + 9 ^{cos (2g) -|- cosh (2/?)). (58) 

Unfortunately one may now derive from m and ()58|1 that the phase shift when the 
two magnons scatter at large L is always as in the BDS chain, and thus as in the gauge 
theory’s near-BMN limit. The exponential corrections disappear in and near the BMN 

limit. 
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4.3 Many Magnons and a Proof of the BDS Formula 

We have seen above that magnons | can arise as bound states of holes o and doubly 
occupied sites The solutions associated to these bound states are known in the con¬ 
densed matter literature as fc—A strings® and were hrst considered by Takahashi® in |27|. 
The explicit solution of the one and two magnon problem allowed us to understand that 
the deviation from the “ideal string” conhguration vanishes exponentially with the chain 
length. In other terms, the string solutions are asymptotic. 

In this chapter we consider the case of solutions with an arbitrary number of magnons. 
We are able to show that, in the asymptotic regime L ^ oo, the scattering of magnons 
associated to the bound states discussed above is described by the BDS ansatz. 

The hnite size corrections may be evaluated, similarly to the one-magnon case, to be 
of the order where jS is the typical strength of the binding sinh/3 ^ 1/g. At weak 

coupling, or in the perturbative regime, these corrections are of order as expected. 

For reasons of simplicity, we are concentrating hrst on magnons with real momentum, 
that is strings containing only one u. In this situation, the momenta Qn appear in complex 
conjugate pairs. Let us choose the labels such that the hrst M momenta have a positive 
imaginary part /9„, while the last M momenta have a negative imaginary part. With the 
experience gained from the one- and two-magnon case we denote 


- 0 = ^ - 1 - y -F , (59) 

qn+M - 0 = Sn 2 + y “ > 0 , 77, = 1, . . . , M . 

where pn will be the magnon momentum, and = sign Pn^^■ If L is large, the left hand 
side of vanishes exponentially for n = 1,..., M and diverges for n = M-|-1,..., 2M. 
Therefore, for L inhnite and for any n = 1,..., M there exist one m, which will be called 
such that 

Un-i/2 = \/2gsin{qn-(t)) , -F i/2 = sin(g„+M - 0) , ( 60 ) 

or, equivalently, 

±i/2 = s„cos(y =Fi0n) • (61) 


In particular, equation (jnj) allows to determine the inverse size of the bound state, /3„, 
in terms of the magnon momentum p^ 

sinh j3n = —^-= — 1 = - (62) 

2\/2g Sn sin ^ 2^/2g\sin^\ 

®In our notation, they should be called q — u strings. 

®In the repulsive case considered by Takahashi, the energy of such a bound state is greater that the 
energy of its constituents, but the wave function is localized in space, so they can still be called bound 
states. 

^°We assume that Pn & (— tt, tt), meaning that the real part oi Qn — (/> ranges from 7r/2 to 3tt/2. It 
is interesting to note that there is no consistent solution with qn — (j) G (—7r/2,7r/2). Such a solution 
would imply a negative energy for the corresponding magnon, which is unphysical. 
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and to find the relation between Un and pn 


Ur 


= V2g 


Pv 


Pr 


COS — cosh (3n = - cot — 
2 ^22 


1 I o 9 • 2 P'^ 

1 + 8g^ sin — , 


( 63 ) 


which is nothing else that the relation m of the BDS Bethe ansatz. The next step is to 
eliminate the fermion momenta from the dnal Lieb-Wn eqnations and replace them 
by the magnon rapidities In order to perform this task, we mnltiply the eqnations 
nnmber n and n + M in dSID, so that the real parts in the exponential mntnally cancel 


gj(p„+20+s„7r)L ^ _ TT Uj + j 

j=l J 


(64) 


Under the condition = (—which is satished dne to onr choice of the twist 
eqnation (IM|) is identical to the BDS Bethe ansatz eqnation CD- The second dnal 
Lieb-Wn eqnation is antomatically satished, while the energy becomes 


M 


\/2 \ M 

E = -> (cosfgn - 0) + cos(gn+M - 0))-y 

^ g 


9 

2 V 2 

9 


(65) 


n=l 

M 


E - Pn 

sm- 


n=l 


cosh0„ - — = - 


9^ 


n=l 


9 ^ 


1 -|- sin^ ^ _ 1 
^ 2 


which is, again, the BDS resnlt ©• 

This proof can be easily extended to the sitnation when the magnon momenta Pn are 
not all real. This may be the case for strings containing more than a single u. We can 
think of snch a string as being composed of several one- magnon strings, each centered 
to a complex momentnm The above eqnations are still valid, nnder the provision 
that Sn is dehned as the sign of the real part of Pn, Sn = sign Re Let ns note that Sn 
is well dehned if Un is hnite. Of conrse, (3n are not real any more bnt they are dehned 
by the hrst eqnality in ©• 


5 Four-Loop Konishi and the Wrapping Problem 

The Hnbbard model is capable of natnrally dealing with the “wrapping problem” [7j. 
The latter is a fnndamental difhcnlty for a long-range spin chain, where one has to 
decide how to interpret the Hamiltonian when the interaction range reaches the size of 
the system^^. 

there are only two-body long-range interactions, as e.g. in the Inozemtsev long-range spin chain 
d, the problem may be circumvented by periodizing the two-body interaction potential. If there are 
also multi-body interactions, as occurs in the long-range spin chains appearing in perturbative gauge 
theory, it is just not clear how to deal with this problem in a natural fashion. See m for a very recent 
discussion of these problems. 
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Let us state the prediction of the Hubbard model for the anomalous dimension of 
the lowest non-trivial state, the Konishi held, with L = 4 and M = 2, to e.g. eight-loop 
order^^. It is easily obtained using e.g. the tool in Appendix iBl 

h^Hubbard = 6-12^2^42c/^-318c/®+4524c/®-63786^^°+783924c/^2_g728086^^^+... (66) 


The four-loop prediction, — 3185 '®, is the hrst order where wrapping occurs. The result 
should be contrasted to the BDS Bethe ansatz, which, when we “illegally” apply it 
beyond wrapping order, yields (again, we used the program described in Appendix lB|l 


h^BPS = 6-12/+425^ 


705 « 6627 




-g 


67287 




16 


359655 

1 ^' 


.12 


7964283 


9^^ + ... (67) 


64 


We can now see explicitly that the perturbative results for the energy differs in the 
two ansatze at 0{g^), i.e. four loop order. The exact result for Konishi is given by a 
rather intricate algebraic curve. Note that the two rapidities ui,U 2 are not related by the 
symmetry Ui = —U 2 - 

Let us likewise contrast the results for the lowest non-BPS state with an odd length, 
namely L = 5, M = 2. The Hubbard model gives 


2 4 115 833 o 6147 44561 ,2 303667 ,4 , , 

h^Hubbard = 4 - 65 ' + 175 " - —g^ + —g^ - -^g^^ + -^g^^ -(68) 


16 


32 


while the BDS ansatz yields 


. . 2 4 115 « 849 „ 6627 

h^BDS — 4 — 65 -I- 175 - —g -\——g - 


5^° + 


53857 
16 ' 


,12 


451507 


32 


+ ... (69) 


In line with expectation this conhrms that the perturbative results for the energy differ 
between Hubbard and BDS at 0{g^), i.e. hve loop order. This is precisely where wrapping 
hrst occurs for a length hve ring. The exact result is again given by an intricate algebraic 
curve. 


6 Conclusions 

The main result of this paper is the identihcation of the long-range BDS spin chain of 
[Zj as an asymptotic approximation to a short-range model of itinerant fermions, the 
Hubbard model. The latter yields a rigorous microscopic dehnition of the former. It 
furthermore provides the Hamiltonian, which was only known, in an “ehective” form, to 
hve-loop order [7] . We have explicitly derived the emergence of this ehective description 
to three-loop order by correcting a previously performed strong-coupling expansion of 
the one-dimensional Hubbard model jlH]. This establishes and proves that the planar 
three-loop dilatation operator of A/” = 4 gauge theory is, in the su(2) sector, generated 
by a twisted Hubbard model. We have also derived the asymptotic Bethe equations of 
the BDS chain from the Lieb-Wu equations of the Hubbard model. 

is interesting that the coefficients seem to be all integer, at least to the order we checked. 
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Our identification allows to resolve the wrapping problem of the BDS chain in a, as 
far as we can currently see, unique fashion. It also gives a rigorous definition of integra- 
bility beyond wrapping order and therefore for a system of finite extent. Recall that the 
notion of “perturbative” integrability implemented in [7j requires, strictly speaking, an 
inhnite system. This renders the BDS ansatz 0 , 0 ,® only asymptotically and thus ap¬ 
proximately valid. The, admittedly more complicated, Lieb-Wu equations (jSni),(ISH),(IS2D 
or (EZD,®,® are the generalization of the BDS equations to strictly hnite systems 
and to arbitrary values of the coupling constant g. Their hrm base is an underlying 
S-matrix satisfying the Yang-Baxter equation ng. What is more, the Hubbard model 
may be included into the rigorous framework of the quantum inverse scattering method. 
In fact, Shastry discovered its R-matrix inoi, and Ramos and Martins inn diagonalized 
the model by algebraic Bethe ansatz. These results therefore also embed the BDS spin 
chain into the systematic inverse scattering formalism. 

We have not been able to hnd the “effective” ansatz 0 , 0 ,®, which signihcantly 
simplihes the nested Lieb-Wu equations at half-hlling up to wrapping terms, in the (vast) 
literature on the Hubbard model uni- This striking simplihcation seems to be a discrete 
and generic generalization of the decoupling phenomenon of the system of thermody¬ 
namic integral equations for the antiferromagnetic ground state energy, as originally 
observed by Lieb and Wu ng. 

Our results strongly indicate that, sadly, wrapping interactions are not able to explain 
the three-loop discrepancies [23 E! between gauge and string theory, as was originally 
hoped for in a proposal in [2|. As discussed in chapter HI the Hubbard model simply does 
not seem to allow for two distinct ways to form the small BMN parameter X' g^Lfi 
Put differently, in the Hubbard model there is no order of limits problem, and wrappings 
just lead to 0{g^^) effects which disappear in the BMN limit. This negative result seems 
to be in agreement with the complementary hndings in m 

Actually, we cannot currently exclude that there might be other, similar (modihed, 
generalized Hubbard?) models which also agree with BDS up to wrapping order, but 
differ from our current proposal in the wrapping terms. However, even if these exist, we 
hnd it hard to believe that they will allow for a new way to form the BMN parameter 
A' at strong coupling. 

These questions should be distinguished from the related, but distinct (since the 
AdS/CFT discrepancy appears at three loops) issue whether the BDS-Hubbard system 
is actually describing the gauge theory’s su(2) dilatation operator at and beyond four- 
loop order. It is of course logically possible that the latter is not asymptotically given by 
the BDS chain at some loop order larger than three. Assuming integrability, we would 
then conclude that BMN scaling should break down at, or beyond, four-loop order, cf. [Z|- 
Then the BMN proposal |H] along with the arguments of ing would be invalid for the 
gauge side. 

It should be clear from the preceding discussion that we are in dire need of a per¬ 
turbative four-loop anomalous dimension computation in J\f = A gauge field theory. Of 
particular importance would be the four-loop dimension of the Konishi field. If it turns 
out to agree with our finding in this paper {—318g^), our attempts to identify the su(2) 
sector of the Af = A dilatation operator with the Hubbard Hamiltonian will, in our opin¬ 
ion, become very plausible. If it disagrees, the search for the correct all-loop dilatation 
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operator will have to be continued. 

Strong additional constraints come from considering the integrable structure of the 
dilatation operator beyond the su(2) sector. The su(2) three-loop dilatation operator j2] 
is naturally embedded in the maximally compact closed sector su(2|3) |3]. The asymp¬ 
totic BDS ansatz may also be lifted in a very natural fashion to this larger sector [HHl I34| . 
Here “natural” means that the ansatz (1) contains BDS as a limit, (2) diagonalizes the 
three-loop dilatation operator in the su(2|3) sector, which is hrmly established |31 E], 
and (3) may be derived from a factorized S-matrix satisfying the Yang-Baxter algebra 
[HSIISIIES]- Actually, the asymptotic BDS ansatz may even be modihed to include non¬ 
compact sectors such as sl(2) [HSl, and lifted to the complete theory [Hll, with symmetry 
psu(2,2|4). Again, the construction seems compelling as it may be shown that (1) the 
Bethe ansatz correctly diagonalizes to three loops twist-two operators jSB] whose dimen¬ 
sions are known form the work of |H7l I38j . (2) it also diagonalizes a twist-three operator 
to two loops which was conhrmed using held theory in [3^]. In fact, it may be proved (3) 
that it diagonalizes to two loops the dilatation operator in the psu(l, 1|2) sector which 
has recently been computed by Zwiebel, using algebraic means in jlH], and (4) for s[(2) 
one may derive the ansatz at two loops directly from the held theory mi. Finally, the en¬ 
tire p3u(2, 2|4) ansatz may again be derived from an S-matrix satisfying the Yang-Baxter 
equation j^. It is important to note that the structure of the S-matrix, as well as, as a 
consequence, the nested asymptotic Bethe ansatz, are nearly completely constrained by 
symmetry IHS] , up to a global scattering “dressing factor” i2niEniin3. This means that 
e.g. the Inozemtsev model |B] is ruled out mi as an all-loop candidate. It also means 
that a possible breakdown of BMN scaling, confer the discussion above, could only be 
caused by the dressing factor, starting at or above four loops. See also j2Hl- Incidentally, 
it would be very interesting to understand whether short-range formulations also exist 
for other (or even all) asymptotically integrable long-range spin chains [3ti|.|28j. 

From the preceding discussion we conclude that it will be crucial to investigate 
whether the twisted Hubbard model may be extended to sectors other than 3u(2), and 
eventually to the full symmetry algebra p3u(2,2|4). A further constraint will be that 
this extension asymptotically yields the Bethe equations of |21]. It would be exciting if 
hnding the proper short-range formulation of the full dilatation operator resulted, when 
restricting to su(2), in a model that also asymptotically generates BDS but differed from 
the specihc Hubbard Hamiltonian we discussed in this paper. At any rate we hnd it likely, 
given the results of this work, that such a short-range formulation of the gauge theory 
dilatation operator exists. It will be interesting to see whether the latter also eliminates 
the length-changing operations which appear in the current long-range formulation as a 
“dynamic” spin chain |3]. 

An intriguing if puzzling aspect of our formulation is that the Hubbard model has 
many more states than the perturbative gauge theory in the 3u(2) sector. For a length 
L operator we have roughly 2^/L cyclic states in the spin chain and in the gauge theory, 
and 4^/L cyclic states in the Hubbard model, cf. section PTTl Is this an artifact of 
the incompleteness, or erroneousness, of our identihcation, or a hrst hint at a rich non- 
perturbative structure of planar A/" = 4 gauge theory? Does it possibly tell us that 
the helds appearing in the Lagrangian of the TV = 4 theory are composites of more 
fundamental degrees of freedom (such as the “electrons” of our model)? A description 
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of the dilatation operator in terms of fermionic and bosonic degrees of freedom akin the 
fermionic variables in the Hubbard model, which works to two loop order, was attempted 
in [321 • Note also that the Hubbard model has a second “hidden” su(2) symmetry jlbj . 
Our twisting procedure actually breaks the symmetry through the boundary conditions. 
Thermodynamically, however, i.e. in the large L limit, the symmetry is still present. The 
mechanism is reminiscent of the considerations of Minahan 123, but the details appear 
to be different. 

Concerning the proposed AdS/CFT duality our result, for the moment, just 
deepens the mystery of the vexing “three-loop discrepancies” of 123,13. The dual string 
theory is classically integrable ^ 2 ] , which leads to a complete solution of classical string 
motions SSI in terms of an algebraic curve HE]. The uncovered integrable structure is 
very similar gn but different i,l3 from the one of the (thermodynamic limit of) gauge 
theory. Much evidence was found that the string theory is also quantum integrable. 
This can be established by a spectroscopic analysis of the spectrum of strings in the 
near-BMN limit |2H|,jlH|, which shows that it may be “phenomenologically” explained 
by factorized scattering 123.123. Again, the integrable structure is similar but, at the 
moment, appears to differ. Some progress has also been made towards deriving quantum 
integrability directly from the string sigma model j49j . 

It would be exciting to hnd a Hubbard-type short-range model which reproduces 
the string theory results. Recently it was demonstrated by Mann and Polchinski isi 
that conformal quantum sigma models can give Bethe equations whose classical limit 
reproduces (in the su(2) sector) the bootstrap equations of |1E|. There is one structural 
feature of their approach which strongly resembles the considerations in this paper: In 
order to be able to treat the su(2) case they need to employ a nested Bethe ansatz, 
which is reminiscent of the Lieb-Wu equations of the Hubbard model. A difference, 
however, is that in their case elementary excitations of the same type are interacting 
with a non-trivial S-matrix, while in our model identical fermions are free. 

Acknowledgments 

D.S. and M.S. thank the Kavli Institute for Theoretical Physics, Santa Barbara, for 
hospitality. D.S. thanks the Albert-Einstein-Institut, Potsdam, for hospitality. We also 
thank Niklas Beisert, David Berenstein, Zvi Bern, Burkhard Eden, Sergey Frolov, Frank 
Gohmann, David Gross, Volodya Kazakov, Ivan Kostov, Juan Maldacena, Nelia Mann, 
Jan Plefka, Joe Polchinski, Radu Roiban and Aureliano Skirzewski for their interest and 
for helpful discussions. This research was supported in part by the National Science 
Foundation under Grant No. PHY99-07949. 


A Effective Spin Hamiltonian: Perturbation Theory 
and Computation Schemes 

A.l Perturbation Theory for Degenerate Systems 

Gonsider a system which is described by a Hamiltonian Hq. Assume that the spectrum 
of Hq is discrete, and that the system is in a stable state with energy E^. In general 
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the subspace Ua corresponding to an eigenvalue has dimension Qa, where ga is the 
degeneracy of the level E^. Let us denote by \ui), ...,\ug^) the vectors spanning Ua- 
What happens if we add a small interaction +W1 In general we have a set of subspaces 
£^i, ...,£n) for which £^i(A) + ...+i^n(A) —Ua when A —0 and dim(£^i + ...+i^n) = da- If A is 
sufficiently small, we may assume that there exists a one-to-one correspondence between 
Ua and W = i^i(A) -|- ... -|-£^n(A). This correspondence is established by a transformation, 
to be found. Let |0) be any state in the Hilbert space generated by E[q. Its projection 
on Ua is formally realized by : 


Pn = 


dz 


2ni Jco z- Ho 


(70) 


where the contour Cq is enclosing only the eigenvalue E^ of Hq. From this discussion we 
conclude that the projector on the subspace W is given by 


1 / dz 

2TTi Jc z — Hq — XV 


(71) 


The contour C encloses the n -|- 1 points Pi(A),..., En{X) (the last n collapse to E^ 
when A —0). Using the identity 


z-Ho-XV z-Ho^^ Ho + XV ^ ^ 

1 1 

z- HZ z- H„-XV’ 


(72) 


we immediately get the expansion 


P= — idz 
2m 


c 


z- Hr 


X;a"(k 


n=0 


z-Ha 


(73) 


Careful use of the generalized Cauchy integral formula leads to the expansion 


P = Po-5^A” S'^^VS'^W...VS’^-+\ (74) 

n=l ki-\-...+knJ^i=n, ki>0 

where one dehnes 

5°^-Po, S^= for k>0. (75) 

Naively one would expect that the correspondence between Ua and IF is realized by the 
projector P, i.e. that any state |'^) G IF can be written as 


W') = PW . 


(76) 


where |0) is some vector in Ua- This would allow us to bring the eigenvalue problem 


P|^) = P|^), 


(77) 
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with 

to the subspace Ua 


H = Ho + XV, 
HP\(j)) = EP\(j)). 


(78) 

(79) 


The disadvantage of this procedure is that EP is not proportional to the identity map. 
Furthermore P does not preserve the norm of the states. The problem of this effective 
overlap has been solved by Lowdin IZQ. One introduces renormalized states (which are 
still states from Ua) 

|0) = (PoPPo)'/"|0), (80) 


and thus one is lead to the introduction of the Ua ^ W correspondence operator T: 


r = pPo(P„pPo )-''^, 


(81) 


where 


(P„PP„)-'/^ - A + E 


n=l 


2n 


n 


p„(p„ - p)p„ 


(82) 


plus an analogous formula for (PqPPq)^^'^. One can then prove that TlT = Pq, so 

{m,m) = m,m), ( 83 ) 


and the transformation preserves the norm. We may now substitute equation dZZI) by an 
effective equation 

(h-P)|0')=O, (84) 

where 

h = T^HT . (85) 

The operator h is the effective Hamiltonian. To hnd it for the Hubbard model at 
half-hlling we put 

L L 

Ho = tU^ oVitoVq, ^ + Xji), Xij = ^ cldr , (86) 

i=l i<j 


One may show that the odd powers disappear from the expansion of h, as one would 
expect from the ’hopping and hopping back’ random walk interpretation: 

h = X^h2 + X^hi + X^ho + .... (87) 


A.2 Computation Schemes 

Performing the computations can be divided into three stages: 
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Stage 1 This stage consists of evaluating the effective Hamiltonian to a given 
order. This is a tedious problem beyond the hrst few orders. One can however use that 
= 0 for the half-hlled Hubbard model, whence in order to get h to n-th order, one 
only needs to evaluate T to (n — l)-th order. Furthermore it can be proved that any 
term of the form 

PoVS^^VS^\..VS^-WPo, ki>l, ( 88 ) 

for odd n vanishes identically. This two observations greatly speed up the calculations. 
A program in FORM (see P) was written to perform this stage of the calculations. 
The result we found up to three-loop reads 

h = +\^{PqVSVPo) + X^PqVSVSVSVPq - \PoVSVPoVSSVPq 

-\P^VSSVP^VSVP^) 

+X^{PoVSVSVSVSVSVPo - IPqVSVSVSVPoVSSVPq 
-^PqVSVSVSSVPqVSVPo - ^PqVSVSSVSVPqVSVPo 
-\P oVSVPoVSVSVSSVPo - \PoVSVP qVSVSSVSVP o (89) 

+lPoVSVPoVSVPoVSSSVPo - \PqVSVPqVSSVSVSV Pq 
+^PoVSVPoVSSVPoVSSVPo - \PqVSSVSVSVPqVSVP o 
-IPoVSSVPqVSVSVSVPo + IPqVSSVPoVSVPqVSSVPo 
+ lPoVSSVPoVSSVPoVSVPo + Xp^VSSSVPqVSVPoVSVPq) , 

It is indeed equivalent to the expansion obtained by Klein and Seitz in [18j . 


Stage 2 This stage consists of substituting (IFFH) into h 2 n as calculated in stage 1. The 
process of substitution can be well visualized by assigning to each Xij an oriented line, 
starting at j and ending in i, see Fig.la. Products of the X operators are represented by 
an oriented set of arrows, with the understanding that the lowest lying arrow corresponds 
to the last operator in the product. A curly bracket around a set of arrows denotes a 
sum over different locations of the arrows. One can interpret these diagrams as virtual 
displacements of spins. It was proved that the perturbation expansion consists only of 
linked diagrams (see HSI for details). Each diagram is multiplied by a suitable factor 
following from the structure of the h 2 n expansion. 


Stage 3 This hnal stage consists of evaluating the diagrams obtained in stage 2. Since 
the diagrams are closed, for each lattice site i the number of arrows starting and ending 
at i is the same. Keeping in mind the dehnition of Xij, and using anti-commutation 
relations for every diagram connecting r lattice sides, we can assign each diagram linear 
combinations of terms of the form 


N{i,,Ti,T 2 )...N(it,T 2 k-,,T 2 t) k<r, 


(90) 


where 

N{i, T, a) = cl^Cia . 


(91) 
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Furthermore one may rewrite each diagram in terms of spin components by means of the 
relations 

St = CpCit, S- = , (92) 

and ^ ^ 

St = - c\^Cii) ~l/2-cJ^Ci^, (93) 

where the last two equalities are only valid when acting on states with a single electron 
per site. This is however our case, after putting the diagrams into the form PI). The 
whole procedure is carried out in FORM. 


a) 


b) 




= Xj,X„ 


FIG. 1. a) To each Xij operator we assign an oriented arrow emerging from site j. b) A 
product of operators is represented by an ordered set of arrows. The lowest lying arrow 
corresponds to the last operator in the product. 


A.3 Twist factors 

Equation m differs from (i2nD by the fact, that each hopping to the right is multiplied 
by a factor of while hopping to the left gets an extra factor of Since the 

perturbation theory consists only of closed diagrams, we conclude that these factors 
cancel at the end. 

This reasoning is generally true for long spin chains. A notable exception is when 
the chain is sufficiently short such that a spin can hop around the ring. This happens 
for example at two loops and L = 4. There are two diagrams corresponding to this 
process. They are related to each other by reversing all arrows in one of them. The 
two diagrams have thus weights differing by factors of opposite signature (e®'^^ = i and 
e-'*^ = c.l Cl) and they therefore cancel each other. Thus putting the twist results 
in eliminating these unwanted demi-wrapping terms. 

B Mathematica Code for the Perturbative Solution of 
the Lieb-Wu Equations 

In this appendix we will collect some Mathematica routines which allow for the immediate 
solution of the BBS equations o,©,© and the Lieb-Wu equations (iHnji . inTTii . in© for 
systems of relatively small lengths L once the one-loop solution is known. The necessary 
input is thus the collection of one-loop Bethe roots {uk}, i.e. the solution of for the 
state in question. The one-loop roots for the first few states may be found in appendix A 
of [ 7 ] . The below routines may therefore be used to check our claims about the agreement 
(below wrapping order) and disagreement (at and beyond wrapping order) between the 
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BDS ansatz and the Lieb-Wu ansatz on various specific states. There is however one 
restriction where the program does not directly apply: There are a number of “singular” 
states nnm with three special unpaired one-loop roots Ui = —^,U 2 = 0,^3 = | which 
require regularization. 

These dehnitions set up the algorithm; 

uu[k_, 1_] := Sum[u[k, slg^s, {s, 0, 21 - 2}]; 

qq[n_, 1_] := Sum[q[n, slg^s, {s, 0, 21 - 2}]; 

\[Phi][L_] := If[EvenQ[L] == True, Pi/( 2L) , 0]; 

x[u_] := u/2(l + SqrtEl - 2 g''2/u~2]) 

BDS[L_, M_, 1_] := 

Table[(x[uu[k, 1] + I/2]/x[uu[k, 1] - I/2])~L + 

Product [(uu[k, 1] - uu[j, 1] + I)/(uu[k, 1] - uu[j, 1] - I), 

{j, 1, M}], {k, 1, M}] 

EBDS[M_, 1_] := Sum[I/x[uu[k, 1] + 1/2] - I/x[uu[k, 1] - 1/2], {k, M}] 

LW1[L_, M_, 1_] :=Table[Exp[I qq[n, 1] L] -Product[ 

(uu[j, 1] - Sqrt[2]g Sin[(qq[n, 1] + \[Phi][L])] - 1/2)/ 

(uu[j, 1]-Sqrt[2]g Sin[(qq[n, 1] + \[Phi][L])] + 1/2),{j, M}] , {n, L}] 

LW2[L_, M_, 1_] := Table[Product[ 

(uu[k, 1] - Sqrt[2]g Sin[qq[n, 1] + \[Phi][L]] + 1/2)/ 

(uu[k, 1] - Sqrt[2]g Sin[qq[n, 1] + \[Phi][L]] - 1/2),{n, L}] + 

Product [(uu[k, 1] - uu[j, 1] + I)/(uu[k, 1] - uu[j, 1] - I), 

{j, 1, M>], {k,l, M}] 

ELW[L_, 1_] := Sqrt[2]/g Sum [Cos [(qq[n, 1] + \[Phi][L])], {n, L}] 

In order to hnd the prediction of the BDS chain for e.g. the state with L = 5 and 
M = 2, where the two one-loop Bethe roots are ui = \ and M 2 = — we then compute, 
to e.g. I = 8 loops, 

Clear[u] ; Clear [q]; L = 5; M = 2; 1 = 8; 
u[l, 0] = 0.5; u[2, 0] = -0.5; 

Do[xxx = Chop[Series[BDS[L, M, 21+1], {g, 0, 11 - 2}]]; 
yyy “ Flatten[Chop[Solve[coefficient[xxx, g, 11 - 2] == 0]]]; 

Do[u[k, 11 - 2] = yyy[[k]] [[2]] , {k, 1, M}] , {11, 3, 2 1 + 1}] ; 

Series[EBDS[M, 21], {g, 0, 21 - 2}] // Chop // Rationalize 
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If we are, on the other hand interested in the correct result of the Hubbard model, 
we compute instead 

Clear[u] ; Clear [q]; L = 5; M = 2; 1 = 8; 
u[l, 0] = 0.5; u[2, 0] = -0.5; 

Do[q[n, 0] = 2 Pi/L(n - 1), {n, 1, L}]; 

Do[xxx = Chop[Series[LW2[L, M, 21+1], {g, 0, 11 - 2}]]; 
yyy “ Flatten[Chop[Solve[coefficient[xxx, g, 11 - 2] == 0]]]; 

Do[u[k, 11-2] = yyy[[k]] [[2]] , {k, 1, M}] ; 

uuu = Chop[Series[LW1[L, M, 21+1], {g, 0, 11-2}]]; 

vvv = Flatten[Chop[Solve[Coefficient[uuu, g, 11 - 2] == 0]]]; 

Do[q[n, 11-2] = vvv[[n]] [[2]] , {n, 1, L}] , {11, 3, 2 1 + 1}]; 

Series[ELW[L, 21+1], {g, 0, 21 - 2}] // Chop // Rationalize 

C Generic Twists 

In this appendix we study all the possible twisted boundary conditions for the Hubbard 
model which are compatible with integrability and the way they affect the Lieb-Wu 
equations. The results are essentially due to Yue and Deguchi who studied the 
twisted boundary conditions associated to a model of two coupled XY models which, 
upon a Jordan-Wigner transformation, is equivalent to the twisted Hubbard model. 
Translating their results in terms of the Hubbard model, we obtain that the twists 
depend on six different constants 


(j)-^ — T X T Ah 
4>i = N + M Cl 


(94) 

(96) 


while the corresponding version of the Lieb-Wu equations is 



H % “ +29 sm((/„ + .)>,) + i/2 ’ 


(96) 



M/j Uj I 

j=l -I 


k = 1, . . . , M 


The energy of the corresponding states is given by 



( 97 ) 
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After the duality transformation, the fermion number becomes L — N -\- 2M, g changes 
sign and > vr — , and ^ 0j,- The dual Lieb-Wu equations are, for generic twist 


^iq-nL _ 


M 

n 


Uj - \/2gsm{qn - 0|) - i/2 


Uj - \/2gsm{qn - 0t) + V2 ’ 


^ Uk - V2gsm{qn - 0t) + 

Uk - V 2 g sin(g„ - 0|) - i/2 




n = l,...,L- N + 2M (98) 

k = 1 ,...,M 


M 

n 

i=i 


Uk — Uj + i 


Uk 


Mi- 


while the energy is 


^ M V2 ^ ^ . 

E = ---> cos(g„ - . 

qz q 

^ ^ n=l 


(99) 


To obtain the BDS ansatz, the following conditions on the twists have to be satished 

( 100 ) 


giL(2*,-,) ^ ^ = Tr{L+l) w 

ZJj 1j 


These are exactly the values we used in m, so we infer that there is no other possibility 
to choose the twists compatible with the BDS ansatz. 


D Further Details on the One-Magnon Problem 

In section 14.11 we discussed how to account for all states of the twisted Hamiltonian 
acting on L — 1 up spins and M = 1 down spin. Recall that in the Hubbard model 
this corresponds to a two-body problem, hence there are states. L of these states are 
bound states, whose dispersion law (I48j) coincides with the one of the magnons in the 
BDS chain. This law turns, using p = 27Tn/L, into the BMN square-root formula 

g"^ E ~ —1 -|- Vl + A' (101) 

if we scale X = 871 “^ g"^ ^ 00 , L ^ 00 while holding A' = A/L^ hxed. As we showed in 
sections o and 01 the scattering of these bound states is as in the near-BMN limit of 
the BDS chain. It is therefore, at third order (A')^, incompatible with the predictions of 
string theory 

One potential way out of this trouble would be to hnd other states in our model which 
scatter as in string theory. A prerequisite is that the coupling constant dependence of 
the dispersion law of such candidate states is again as in (HOD), with, possibly, a different 
constant part. In particular, among the real solutions we identihed in .section ITTI there 
were states of exact dimension A = L —1 which resembled “BPS states”. Let us therefore 
work out the dispersion law of the nearby “near-BPS” states. This requires studying the 
solutions of (EH) for small q = ir/Ln and large g. Expressing (IKK]) through the BMN 
coupling A', we hnd 

sin (6) = _ tan(6L) . (102) 

vA'n 
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Since this equation should hold as A' —0, we recover the L — 1 mode numbers m (jSl: 


b 


71 

— m + 



(103) 


Now, 6 bm should be at most of order 0(1), i.e. it should not be too large so as to move 
out of the branch of tan{bmL) dehned by (ll(Kf|l . and should tend to zero if y/}/ u —>• 0. 
This yields from (uni 5bm — y/^n sin(^m). Substitution into the expression for the 
energy of the real solutions 


E 


gives for the dimension A 



2V2 

9 


sin (g) cos(6) 


A 


, TT 


L — 1 + y/y n cos(— m) 

1j 


^ sin^ + 0{1/L) . 


(104) 


(105) 


We see that we generically lift the L — 1 degenerate “BPS-states” with a term non- 
analytic in A'. If we concentrate on mode numbers close to m ~ L/2 we can suppress 
the non-analytic y/)^ term. The next term is then analytic in A', but subleading in 1/L. 
It is interesting to note that there is a possibility to reproduce a BMN-like dispersion 
relation, by choosing m such that 

cos(y m) =-^V^n + (P(l/L) , (106) 

1j z 


so that the conformal dimension would be analytic in A' up to terms of order 1/L 

A = L - 1 + ^A'n^ + 0{l/L) . (107) 

However, such a choice for m is not continuous in A' and cannot be sensibly interpreted 
in terms of BMN states. The “BPS-states” we found are thus very different from the 
usual ones, and the BMN states may not be expected to hide among the continuum of 
real solutions. 

Finally note that any one of the L bound states of section IT?T] can disappear^^ if g 
is very close to L. One may show that in this case a further real solution with mode 
number m = 0, which generically does not correspond to a solution of (ESI) , appears. 
Unfortunately this deconhnement phenomenon is also not suitable for hnding the BMN 
states of string theory [Hj, as we are then not allowed to make the parameter in (j1()2j) 
arbitrarily small. 


E Alternative Proof of the BDS Equations 

In this appendix we give an alternative proof to the BDS ansatz, using the original Lieb- 
Wu equation, with a macroscopic number of fermionic excitations. As in the original 

^^This is the so-called “redistribution phenomenon” m and is responsible for rendering all the fermion 
momenta real in the extreme g limit, g ^ L, which corresponds to the free fermion limit. 
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paper HSI, we suppose that the fermion momenta are all real and they form, in the 
continuous limit, a continuous density. This proof is less effective than the one which 
starts from the dual Bethe ansatz, in the sense that the finite size corrections are not 
under control, and the effect of the boundary conditions (twist) is lost. However, it is 
interesting to see that the BDS equations are already contained in the integral equations 
of Lieb and Wu ini. 

At half-hlling, the Lieb-Wu equations can be written in the logarithmic form as 


27m 2 

qn = (p-\ -7-y 2^ arctan 


L 


2 arctan 


i=i 

1 


2{uj - V^gsmqn) 


M 


, n = 1,..., L 


n=l 


2{uk - V^gsinqn) 


271171 + 2 arctan 


7=1 


2{uk - Uj) ■ 


(108) 

(109) 


The choice of the branch of the logarithm in (HHHi) is made by continuity, such that at 
g = 0 there is exactly one electron per level. For simplicity, we have remove the tilde on 
the variables and shifted them by 0. Taking the derivative of the first equation with 
respect to q and dehning the density p{q) = {dn/dq)/L we obtain an equation for the 
density 


27ip{q) 



2 \/2gcosq 

4{uj — V2gsmqy + 1 


( 110 ) 


Our purpose is to study the case of a finite (arbitrary) number of magnons, so we do 
not introduce a density for the magnons. instead, we evaluate the left hand side of the 
second equation Lieb-Wu equation (ITO 


I{uj) =2L dq p{q) arctan 


2{uj 


1 

\/2g sin q) 


( 111 ) 


The second term in the density does not contribute to the integral I{u). To compute 
the integral /(u), we hrst take its derivative with respect to u, so that the cuts of the 
integrand disappear 


d L 

—liu) =- 

du 2771 


dq 


u 


+ i/2 — \/2gsm.q u — i/2 — y/2gsmq 


( 112 ) 


The integral over q can be traded to a contour integral by a change of variable z = 
\/2g sin q 



L f idz / I 1 \ 

271* fc - 2g^ u--z) 


(113) 


where C is the contour encircling the interval [—\/2g, V^g] clockwise. The contour C 
cannot be shrunk to zero because of the obstruction created by the square root in the 
integrand. The integral vanishes on the contour at inhnity, so we can deform the contour 
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C into two contours C+ and C- which encircle the points and u counterclockwise. 
We obtain 


^I{u) = -iL ( , ^ ^ 

du I a/m+ - 2g^ ^u- - 2g^ 


d x(m+) 
^ du x(u-) 


(114) 


The constant of integration can be easily seen to be zero, since I{oo) = 0. Finally, the 
second Lieb-Wu equation (fTO takes the form 


/ x+{uk) \ _ A -Uj + i 

\x-{uk)J fJiUk-Uj-i' 


(115) 


The magnon energy can be computed by the same means. In this case, only the second 
term in the density (Hini) contributes 


E 


[ dq p{q) cos(g) 


^ 2y/2gcos‘^q 

^9 ^ 4:{uj - V2g sin qy + 1 


(116) 


Again, the integral can be converted into a contour integral around the same contour C 
which encircles the cut [—\/2g, clockwise 


E 


I ^ 

Jc 27r* {z-uj){z- uj) ■ 


(117) 


As such, the integral does not vanish on the contour at inhnity, but we can freely add 
to it a term which is regular across the cut and which removes the contribution from 
inhnity 


iz-up{z-u~) 

I dz x-\z) -,f^( ^ _ L_\ 


(118) 


Of course, the reader recognizes (HinD and (HIHD as the equations of BDS ansatz. 
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